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SUMMARY 


Results of an experimental investigation of propeller whirl flutter are 
presented for a model consisting of an isolated, rigid system of propeller and 
simulated power plant mounted with flexibility in pitch and yaw on a rigid sting. 

A range of propeller blade angles, restraint stiffnesses, and restraint damping 
coefficients was investigated for a system symmetrical in pitch and yaw with a 
windmilling propeller. Measurements of the static-stability derivatives were 
also made by using a simple balance and were compared with two sets of theoretical 
derivatives. Whirl-flutter calculations were made with the theoretical and meas- 
ured derivatives. Some limited results were obtained for the whirl flutter of 
the model mounted on a cantilever semispan wing. 

The measured whirl-flutter speeds and frequencies of the isolated model were 
in very good agreement with those predicted by calculations in which measured 
derivatives and viscous damping were used. This agreement was better than that 
obtained by using structural damping. Predicted whirl- flutter speeds for the 
isolated model were lower when theoretical stability derivatives were used than 
when measured derivatives were used. The theoretical and experimental static- 
stability derivatives exhibited the same trends, but in certain instances dif- 
fered appreciably in magnitude. There was little effect of the semispan wing on 
the measured whirl-flutter boundary for the one configuration considered. 


INTRODUCTION 


It is well known that a combination of propeller and power plant can develop 
a dynamic instability, generally termed propeller whirl flutter or autoprecession, 
in which the hub wobbles or executes a whirling motion. In several analytical 
investigations (refs. 1 to 4) the flutter boundaries of an isolated nacelle have 
been investigated by developing the equations of motion with calculated stability 
derivatives for the aerodynamics of the propeller, and by making parametric trend 
studies from these equations. However, there have been no comparable experimental 
data available for verification of these analytical studies. For example, in 


reference 2 a comparison of the theoretical results with the flutter data for a 
four-engine model showed comparable trends but poor quantitative agreement, 
perhaps as a result of the effects of the complete model on the whirl-flutter 
boundaries. Therefore, the present investigation was undertaken with two 
objectives: 

1. To obtain experimental whirl-flutter boundaries on a simple system that 
is readily amenable to analysis, in order to evaluate existing propeller-whirl 
theories . 

2. To measure, on the same system, some of the propeller aerodynamic deriva- 
tives required to predict whirl flutter, in order to evaluate methods for cal- 
culating the stability derivatives of the propeller. 

Also, the effects of a wing on the flutter boundary are indicated by a com- 
parison with some limited results obtained with the model mounted on a cantilever 
semi span wing. 

Measurements of the static aerodynamic stability derivatives of the wind- 
milling propeller were made by using a simple balance and were compared with the 
theory of reference 3 and the theory of reference 5 as applied in reference 2. 
Whirl-flutter speeds and frequencies were measured for a combination of a wind- 
milling propeller and simulated power plant which was mounted on a sting by means 
of a spring-restrained gimbal. A range of propeller blade angles, restraint 
stiffnesses, and restraint damping coefficients was investigated for a system 
symmetrical in pitch and yaw. Whirl-flutter calculations made with both sets of 
theoretical derivatives (theories of refs. 3 and 5) and measured derivatives 
were compared with experiment. 


SYMBOLS 


A 

b 


propeller-blade aspect ratio 
propeller-blade section chord, ft 


-'in 


pitching-moment coefficient, F L tchi^ _ momen t 


Cm 


total pitching-moment coefficient, C m - — Cg 

2R 


'n 


y awi ng -mome nt coefficient, 


Yawing moment 
qSD 


^n 


total yawing-moment coefficient, C n + — Cy 

2R 


°n 


total yawing-moment coefficient including effect of yaw strain-gage 
beam deflection (see eq. (A7)) 
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Side force 

qS 


side-force coefficient. 


vertical- force coefficient. 


Vertical force 

qS 


section lift-curve slope of propeller blade 


propeller diameter, ft 
vibration frequency, cps 
structural damping coefficient 


average structural damping coefficient, 


S 9 + Pq 

2 P 


moment-of -inertia ratio. 


jtlx 

lY 


mass moment of inertia of propeller about axis of rotation. 


slug-ft 


2 


mass moment of inertia of propeller-engine system about pitch axis, 
slug-ft^ 

propeller advance ratio, — 

’ nD 


rotational stiffness of engine-propeller system, ft-lb/radian 
reduced frequency of engine-propeller system, — 


distance from plane of propeller-blade quarter-chord points at three- 


quarter radius to gimbal axes, 2 0 + 


3 0.75R 


sln e 0.75B’ ft 


distance from plane of propeller-blade midchord line to gimbal axes, 
ft 

free-stream Mach number 
propeller rotational speed, rps 
free-stream dynamic pressure, lb/sq ft 
propeller radius, ft 
propeller disk area, jtR^, sq ft 



time, sec 

free-stream velocity, ft/sec 

distance along propeller blade, fraction of R 

propeller geometric blade angle measured from plane of propeller 
rotation, deg 

propeller blade angle for zero lift measured from plane of propeller 
rotation, deg 

lag angle of oscillatory aerodynamic forces, deg 

viscous-damping coefficient, ratio of viscous damping to critical 
damping 


average viscous-damping coefficient. 


^0 + Po 


2 p 


pitch angle of propeller- shaft axis 


density-inertia ratio, 


frequency ratio, — 

CD 

air density, slugs/cu ft 

reference air density, 0.0022 slug/cu ft 

8b 0 75R 

propeller solidity at three-quarter radius, — ' i; (for four-blade 

propeller) 

nondimensional distance traveled in propeller radii, ~ 


yaw .angle of propeller-shaft axis 
propeller rotational speed, radians/sec 
frequency, radians/sec 


CD q + CD, 

average wind-off frequency, — — % radians/sec 
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Subscripts: 


0.75R three-quarter propeller radius 

f condition at whirl flutter 

i = 1, 2,3 wing vibration modes 
q differentiation with respect to 0 1 

differentiation with respect to \j/ ' 

pitch direction 
yaw direction 

Primes denote differentiation with respect to t- 

Partial derivatives are denoted by double subscripts; for example: 
c>CL 


r 

9 


m 


" W’ 


■ / m r 


m 


dq 


, and so forth. 


APPARATUS 


Tunnel 

The Langley transonic dynamics tunnel is a slotted- throat, variable- 
pressure, single-return wind tunnel having a test section 1 6 feet square (with 
cropped corners). It is capable of operation at Mach numbers up to 1.2 and at 
stagnation pressures from near vacuum to slightly above atmospheric. Either air 
or Freon-12 can be used as a test medium. Large windows are provided for close, 
unobstructed viewing of the model. 


Models 

The isolated model consisted of a combination of simulated engine mass and 
propeller attached to a steel mounting beam through a gimbal with pitch and yaw 
freedoms. The _ engine-propeller system was elastically restrained for pitch and 
yaw movements with a steel spring. The steel mounting beam was bolted to a sting 
on the wind-tunnel center line. A balsa nacelle cover, similar to one type of 
present-day turboprop nacelle, covered the gimbal, propeller shaft, and engine. 
This nacelle was faired at its after end into an aluminum cylinder which covered 
the end of the sting. Figure 1 is a photograph of the model in the tunnel. Two 
propeller-shaft lengths (shown in fig. 2) were used. Details of the method of 
spring restraint and the gimbal arrangement, which differed for the models used 
in the whirl and derivative tests, are given subsequently. 
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The whirl-flutter model was also mounted on a cantilever semispan wing as 
shown in figure 3- 

Propeller .- The four -blade aluminum propeller used in these tests was a 
1 /8 -scale model of that treated in reference 2 and was typical of the type used 
on turboprop aircraft; however, the propeller did not have blade-root cuffs. 

Figure 4 gives the blade width and twist distributions. Figure 5 presents the 
variation of windmilling-propeller advance ratio with blade angle at the three- 
quarter- radius station as found in this test. The propeller blade frequencies, 
with the hub rigidly clamped, were as follows: First bending, 172 cps; first 

chordwise bending, 296 cps; second bending, 567 cps; first torsion, 1,105 cps. 
There was little change in the frequency of the first mode with the propeller 
mounted on the model. 

Isolated whirl model .- Figure 6 is a photograph of the whirl model with the 
cover removed; the model is shown schematically in figure 7 ( a )- The gimbal 
allowed the propeller-engine system freedom to pitch and yaw. The spring 
restraint was a hardened piece of l/ 4 -inch-diameter drill rod, screwed into the 
rear of the propeller shaft housing on the center line of the shaft axis. The 
opposite (rearward) end of the spring passed through a hole in the steel clip 
that could be screwed to the mounting beam at several locations to provide varia- 
tions in stiffness. This arrangement provided equal stiffness in all directions. 

A variable- speed motor with an eccentrically mounted weight on its shaft formed 
part of the engine mass and served as a shaker device. 

The measured physical properties of the model are given in table I. Measured 
values of spring stiffness agreed with values calculated from the measured moments 
of inertia and frequencies to within 5 percent. The pitch and yaw natural fre- 
quencies (propeller not rotating) ranged from 7-6 to 11. 5 cps. The amount of 
damping in the system was adjusted by inserting sponge rubber between the gimbal 
rings (fig. 6 ). Figure 8 shows four typical samples of the data reduced from 
decay records that were obtained by locking out one degree of girafcal freedom, 
deflecting the propeller shaft, and releasing it. In two of the cases presented 
the damping is reasonably independent of amplitude (constant slope); in the other 
cases there is an increase in damping with increasing amplitude. 

For the whirl tests, the sting was restrained by cables attached to its 
upstream end. (See figs. 1 and 6 .) The sting stiffness was about one hundred 
times the stiffness of the propeller-engine system. The sting pitch frequency 
with the model and cables in place was 23 cps; the yaw frequency was 20 cps. 

Derivative model . - The isolated whirl- flutter model was modified in several 
ways for use in the derivative measurements. One direction of gimbal freedom 
was locked out so that the gimbal became essentially a single pivot and allowed 
motions in one plane only as shown in figure 7 (b)- Also, a flat spring was used 
as a strain-gage beam bolted at its forward' end to the propeller housing and 
clamped at its rearward end to a flexible steel clip bolted to the mounting beam. 
The mounting beam could be bolted to the sting in such a manner as to allow 
either pitch or yaw freedom in the tunnel (fig. 7 (b)). 
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It was impractical to restrain the sting with cables during the derivative 
tests, as the end of the sting translated somewhat when the sting angle was 
varied to determine the static derivatives. The sting stiffness was about twenty- 
five times the spring stiffness, and its frequency with the model in place was 
about 14 cps. 

Wing-mounted model . - The propeller and gimbal arrangement for the whirl- 
flutter model was mounted on a cantilever semispan wing. The general arrangement 
and the wing modal frequencies and node lines are indicated in figure 9* The 
wing was adapted from that of the four-engine model of reference 6. It consisted 
of a built-up aluminum spar for stiffness and balsa pods for airfoil contour. 


Instrumentation 

Tunnel static and stagnation pressures were measured with manometers and 
were recorded in digital form, with stagnation temperature, on manual command at 
each data point by means of the tunnel automatic-readout system. 

The angular accelerations about the gimbal axes of the propeller whirl model 
were determined by linear accelerometers mounted on the simulated engine mass as 
shown in figure 6, and were recorded on a direct-writing oscillograph recorder. 
Propeller rotational speed was measured by a magnetic pickup mounted at the pro- 
peller shaft and driving an electronic counter. 

In the stability-derivative tests, the moment about the gimbal axis was 
sensed by strain gages on the gimbal restraining spring. The signal from the 
strain gage was read out on a fixed-gain, digital voltmeter, and printed out in 

digital form. Propeller angle of attack was measured by a ig accelerometer 

mounted on the simulated engine mass (fig. 7(b)). The mean value of the accel- 
erometer signal was recorded in a manner similar to that for the moment signal. 
The sting angle was also measured by a calibrated counter on the actuator mechan- 
ism and was recorded manually. 


TESTS 


Whirl -Flutter Tests 

The linear theory of the stability of the whirl mode (refs. 2 and 3) indi- 
cates that there is a critical speed at which a disturbance of the mode will 
produce constant -amplitude oscillations; at lower speeds the mode would be 
stable and at higher speeds the mode would be unstable. In the tests of the 
whirl- flutter model, the tunnel was brought up to a low speed, the model was 
disturbed by the shaker, and the stability of the response was observed visually 
and on the oscillograph record. Airstream velocity was then increased by a 
small increment, and the model was again disturbed. This procedure was continued 
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until a sustained, approximately constant amplitude oscillation was produced. 

The tunnel conditions were then recorded as the conditions for flutter. 

In some cases, the stability of the whirl mode was related to velocity as 
indicated by linear theory; however, in many of the tests there was no such 
sharply defined flutter speed. In the latter cases the amplitude of flutter was 
limited for a given velocity and increased with an increase in velocity, indi- 
cating a nonlinear effect. These configurations corresponded to those having a 
damping that increased with amplitude (fig. 8 ). For these configurations, veloc- 
ity was increased until the amplitude of oscillation was about 3° or 4°, and- the 
tunnel conditions were then recorded as the conditions for flutter. The pertinent 
damping coefficient used in this report for these cases was assumed to be that 
derived from the slope of the decay records (fig. 8 ) at the amplitude of flutter. 
Damping records were taken before and after a series of runs for each stiffness 
or damping level. 

All the whirl- flutter tests were made in air at near atmospheric conditions 
and low Mach numbers (M ^ 0.2 5 ). 


Derivative Tests 


The moment about the gimbal axis was measured over a wide range of pitch 
angle 0 by varying the angle of the sting. Two distances between the gimbal 


axis and propeller plane were used 


{¥ - °- 51 ‘ 6 


and 



in order to 


separate the force and moment contributions to the total moment. The method of 
obtaining the stability derivatives from the slopes of the total -moment- 
coefficient curves is given in appendix A. 


At low Mach numbers (M % 0.26) the moment about the gimbal axis was measured 
with J3q = 25 °, 35 °, 46°, 52°, and 58 °, with both shaft lengths, and with the 

gimbal axis oriented for pitching-moment measurements and for yawing-moment meas- 
urements. The nominal value of dynamic pressure was either 50 or 100 pounds per 
square foot. In addition, measurements were made for 0O.75R = 58°, and 

— = 0.346 with the gimbal oriented for pitch measurements at Mach numbers of 
R 

O. 9 O, O. 7 O, and 0.48 and oriented for yaw measurements at Mach numbers of 0.90 
and 0.70. A nominal value of dynamic pressure of 100 pounds per square foot was 
maintained at the higher Mach numbers by varying the static pressure within the 
tunnel . 


During some of the derivative tests a mild oscillation of limited amplitude 
was encountered. In some cases it occurred, only at the higher angles of attack, 
and in others at the higher values of dynamic pressure. These oscillations had 
no apparent effect on the measurement of the static derivatives, since such meas- 
urements depend primarily on the mean values of the measured quantities. Attempts 
were made to measure the damping derivative Cmq by obtaining decay records after 

a shaker disturbance at several airstream velocities for a system with a weak 
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pitch spring and no yaw freedom. However, again limited-amplitude oscillations 
were obtained, and thus Cmq could not be determined in a simple manner by using 

the equation for a single degree of freedom. (A method of estimating Cmq from 

whirl-flutter speeds is described subsequently. ) The ratio of pitch to yaw stiff 
ness for this condition (sting tied down with cables) was approximately 1 to 100, 
whereas for the static derivative tests (sting not tied down with cables) the 
stiffness ratio was approximately 1 to 25- The linear two- degree- of- freedom 
analysis of reference 2 for a pivoting simple system does not predict an insta- 
bility at the speeds of the derivative tests for the stiffnesses involved, and 
the explanation of these limited-amplitude oscillations is unknown at the present 
time. 


RESULTS AND DISCUSSION 


Presentation of Results 

The results are presented and discussed in three parts. First, the theoret- 
ical and experimental stability derivatives are considered (figs. 10 to 13) 
because they are required as an input for the whirl- flutter calculations. Second 
the measured and calculated whirl-flutter boundaries and frequencies for the 
isolated nacelle are considered (figs. 14 to 20 and table II ). Third, some wing 
effects are discussed by comparing the measured flutter boundaries for the iso- 
lated nacelle with limited results for the wing-mounted nacelle (fig. 21 and 
table II). 


Propeller Stability Derivatives 

The experimentally determined static-stability derivatives (M ^ 0.2 6) are 
compared in figure 10 with derivatives calculated by the method of reference 5 
(incorporating the lag terms of ref. 2), and with the strip theory of refer- 
ence 3. The derivatives of figure 10 are derived from the faired curves of fig- 
ure 11 by means of equations (a4), (A 5), (A9), and (AlO) of appendix A. The 
theoretical values of C^q (fig- 10(c)) were computed from Cm^ tan 5, which 

can be derived in a manner similar to the derivation of in reference 2.. 

As previously discussed, no experimental values of C m ^ were determined. 

Application of theory .- In the application of the theory of reference 5 to 
the windmilling propeller, the effective local helix angle was assumed to be 
given by V/ftRx. That is, momentum losses and interactions between the blades 
were assumed to be negligible. A section lift-curve slope of 2n, fineness ratio 
of 6.0 for the spinner-nacelle combination, and a ratio of spinner radius to 
propeller radius of 0.l8l were used in evaluating the integrals over the portion 
of the blade having lifting sections (outer three-quarters of the radius). 

The derivatives calculated by the theory of reference 3 for windmilling 
propellers are somewhat simpler to evaluate than those of the theory of 
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the 


reference 5* 


c v 


Since these derivatives are directly proportional to 

^ (developed for wings) was used as sug- 


aspect-ratio correction factor 

gested in reference 3. Evaluating A over the outer tlqree- quarters of the radius 
yields a value of 3.V7 for A and a correction factor of 0.635* The integrals 
of this method (eqs. (B9) of ref. 3) were also evaluated over the outer three- 
quarters of the radius. 


Comparison of theory and experiment . - In general, the trends of the theo- 
retical derivatives correspond to the experimental trends, but the magnitudes 
differ considerably in certain cases. The values of (fig* 10 (a)) predicted 

by the method of reference 5 are about 25 percent larger than the experimental 
ones, while the values predicted by the method of reference 3 are in agreement 
at low blade angles but are smaller by 15 percent at the higher blade angles. 

In references 5 and 7 the data cited were within ±10 percent of values of 

calculated by the method of reference 5. The difference in the present case 
might be attributed to experimental inaccuracy (estimated to be ±10 percent), to 
the value 2« used for section lift-curve slope (section data were not available 
for this propeller), or to the peculiar geometry of the inboard portion of the 
propeller. 


The values of C m q (fig* 10 (a)) computed from the two theories have the 

same trend, but the values predicted by the theory of reference 3 are somewhat 
smaller. Both theories give results for Cm^ that are in good- to -fair agree- 
ment with experimental results over the blade-angle range (fig. 10(b)), but are 
2 to 20 percent higher. The experimental accuracy of this derivative is esti- 
mated to be approximately ±10 percent also. 


The theoretical values of and (fig* 10 (c)) are lower than exper- 

iment and are in fair-to-poor agreement. The estimated experimental accuracy is 
about ±10 percent for 0^ and ±25 percent for C^. The theoretical expres- 
sions for Cjdq and Cz^ (Cmg = Cm^ tan 5; Cz^ = -CZq tan &), as applied in ref- 
erences 2 and 3, are approximations and would not be expected to yield accurate 
results. 


A further indication of the error in the unsteady- flow lift -lag phase 
angle B, as obtained in references 2 and 3, can be seen in figure 12 by com- 
paring the phase angles derived from the experimental derivatives with those 
derived from the theory of reference 2. Although the phase angles derived from 
the experimental force and moment relations should agree, the difference between 
the two experimental curves is within the estimated accuracies of the derivatives. 
In view of the previously cited accuracy of Cz^, ±25 percent, the phase angles 

derived from the moment relations are considered more nearly accurate. 


The flutter calculations of this report are based on the total moment 
derivatives: 
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c 


ze 


~ ~ C nQ “ % - % 

These experimental total moment coefficients are shown in figure 11 . The esti- 
mated accuracies for these quantities are ±5 percent for and ±10 percent 

for C^. 

Mach number effects .- The effect of Mach number on the total moment coeffi- 
cients and for P Q ^ ^ = 58° is shown in figure 13. Little effect 

is indicated for C^q, but appears to decrease with increasing Mach number. 

In contrast, theory would predict an increase in these derivatives given approxi- 
mately by the Prandtl-Glauert correction factor ^ at the lower Mach 

numbers . 


Propeller Whirl Flutter 

Whirl-flutter speeds .- In figure 14 the experimental whirl-flutter data for 
1 0 /R = 0.3^6 are compared with theoretical boundaries for five blade angles. 

The data are presented in the form of damping required for stability as a func- 
tion of reduced velocity V/Hs, the stable region being above the boundary. Since 
the air density varied somewhat from one test point to another (table II), the 
stability boundaries in figure 14 - have been adjusted to represent a reference 
density value p Q of 0.0022 slug/cu ft. As shown in appendix B, the appropriate 

density correction factor is P 0 /p* Also, an average value of damping in the 

pitch and yaw directions is taken to be the representative damping, as refer- 
ence- 2 indicates no effect of the ratio of pitch damping to yaw damping for a 
symmetrical system. Thus, the damping parameter plotted is 

- _ fajjjr Po 

g " 2 P 

or 

z _ + £\|r Po 

^ 2 p 


As considerable scatter exists in the damping measurements, a faired curve was 
obtained by fitting all the data for 2 0 /R = 0*3^6 in the least-squares 
sense to determine the coefficients of the following expression: 
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(123.2 - l. 59^00. 75R 




This particular form can be justified as follows: The approximate whirl equa- 

tion (eq. (26) of ref. 2) indicates that the theoretical whirl boundary may be 
approximated very closely by a parabola in V/Roj. A plot of the coefficients of 
this equation (eq. (26) of ref. 2) shows that they vary in a nearly linear man- 
ner with blade angle. 


Figure 15 presents a comparison of the experimental and theoretical whirl- 
flutter boundaries for i 0 /R - 0.691* Since all the measurements were made at 

approximately the same value of damping, 2£ = 0.0059 is taken as an average 
value and the data are presented in the form of reduced flutter speed v/Rcd as 
a function of blade angle. The stable region is below the curves. Observation 
of the data suggested the use of a parabola to fit the data. Therefore, the 
faired curve in this case was obtained by fitting the data in the least-squares 
sense to determine the coefficients of the parabola: 

^ = (3,054 - 6.6ip 0 . T5R + 0.5243pg. 75R )lO-3 

Theoretical flutter calculations would predict a curve of more complicated shape, 
but this approximation was considered satisfactory for the limited amount of 
data available. 

Theoretical whirl boundaries were calculated by using equations (b 6) and 
(B8) . Since the nature of the damping in the system is unknown, calculations 
were made with the assumption of either viscous or structural damping. In these 
calculations the measured static- stability derivatives and C m q calculated by 

the method of reference 5 were used. For this system with l 0 /R = 0.346 
(fig. 14) and l Q / R = O.691 (fig. 15) the use of viscous damping gives better 

agreement with experiment in that the errors axe generally smaller and the shapes 
of the experimental and theoretical curves are more nearly the same. The agree- 
ment between the theoretical and experimental curves with viscous damping is 
excellent for l 0 /R = 0.346, the largest error in reduced flutter speed being 

6 percent. The error in speed for Z 0 /R = O.69I is generally less than 
10 percent. 

Although the use of viscous damping in the calculations provides better 
agreement with measured whirl speeds in this case, such a result may not be true 
in general, particularly for aircraft structures. 

Figures 16 and 17 show a comparison of experimental and theoretical 
(eqs. (B6) and (b 8)) whirl boundaries for the measured static- stability deriva- 
tives (with the Cinq of ref. 5j, the derivatives of reference 5, and the deriva- 
tives of reference 3* Viscous damping is assumed. It is seen that the use of 
measured derivatives provides the best agreement. The calculated derivatives of 
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reference 5 provide better agreement than those of reference 3’ The use of 
theoretical derivatives is conservative in that lower whirl speeds are predicted 
than were measured. 

Figure 18 gives a comparison of the theoretical whirl speeds for 
Z Q / R = 0.346 and Zo/R = O. 69 I for three blade angles. It is seen that 

increasing Zo/R leads to a more stable system (higher flutter speed) , as was 
indicated by references 2 to it-. 

Cmq required to predict whirl- flutter speeds.- Since no experimental values 

were available, it was of interest to determine the C m q required to predict 

the measured whirl speeds. This was done by using the measured system properties 
of inertia and viscous damping coefficient with the measured whirl speeds and 
measured static derivatives in equations (b 6) and (B8) to solve for C m q. 

Ideally, this might also be considered the Cmq inferred from the whirl-flutter 

speeds, but such a procedure should be considered to yield only an estimate in 
this case. The values of C m computed for two values of damping when 

SI 

Zo/R = 0.346 and for the single value of damping when Zo/R = O. 69 I are com- 
pared with theoretical values in figure 19* The inferred values of Cmq for the 
three cases considered are in good agreement except for the data for Zo/R = O. 69 I 

at the higher blade angles. The deviation of these points can possibly be attrib- 
uted to the sensitivity of the whirl speeds to restraint damping, as indicated 
by the small slope for the higher blade angles at Z 0 /R = O. 69 I in figure 18 . 

The inferred values of C m q (fig. 19) exhibit the same trends as the theoretical 

values and are in good agreement with the values of reference 5 over most of the 
blade-angle range. 

Whirl frequencies .- Figure 20 gives a comparison of experimental and theoret- 
ical whirl frequency ratios for Zo/R = 0.346 and for Zo/R = O. 691 . The theory 
is shown to three degrees of approximation: first, including first-order aero- 

dynamic terms (eq. (B6)); second, neglecting aerodynamic terms ^C^q = 0 in 

eqs. (B7)); and third, neglecting aerodynamic terms and assuming the angular- 

momentum ratio ^2^ « 1 (eq. (BIO)). It is seen that the calculations in which 
Iyco 

C^q is retained provide the best overall agreement, the largest error being 
8 percent. 

Wing effects .- Some tests were made in which the nacelle- engine -propeller 
model treated in this paper was mounted on a semispan wing cantilevered from the 
tunnel wall as previously discussed in the section entitled "Models." Whirl- 
flutter speeds for two levels of damping are presented in figure 21, and are 
compared with results from the tests of the isolated model. This comparison 
indicates that the wing had little effect on the whirl speeds. Whirl fre- 
quencies ranged from 5*7 C P S to 8.9 cps. 
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The question of the effects of the wing is a broad problem, and therefore 
the small effect illustrated in figure 21 may not be typical of the results for 
other configurations. It should also be noted that in the experimental, results 
presented in reference 2 for a four-engine model, the presence of a wing and 
complete model have large effects (up to 100 percent) . 


CONCLUSIONS 


A wind-tunnel investigation of propeller whirl flutter was made with an 
isolated combination of a windmilling propeller and simulated power plant over 
a range of propeller blade angles, restraint stiffnesses, and restraint damping 
coefficients for a system symmetrical in pitch and yaw. Measurements of the 
static- stability derivatives of the windmilling propeller were also made with a 
simple balance. Limited whirl- flutter results were obtained for the model 
mounted on a cantilever wing. The following conclusions are indicated: 

1. The measured whirl- flutter speeds and frequencies for the isolated model 
were in very good agreement with those predicted by calculations in which meas- 
ured static derivatives and viscous damping were used. This agreement was 
better than that obtained using structural damping. 

2. Predicted whirl- flutter speeds for the isolated model were lower when 
theoretical stability derivatives were used than when measured static derivatives 
were used. 

3. Theoretical static- stability derivatives obtained by two methods exhib- 
ited the same trends as experimental derivatives, but differed appreciably in 
magnitude in certain instances. Limited measurements at one blade angle indi- 
cated that Mach number had little effect. 

4 . Little effect of the wing on the measured whirl-flutter boundaries was 
indicated by the limited data available for the one configuration considered, 
although previous investigations have shown larger effects of the wing. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 8, 1 963 . 



APPENDIX A 


DATA REDUCTION FOR STATIC-DERIVATIVE MEASUREMENTS 


Determination of C mg and 


The static equation governing the pitch deflection of the engine-propeller 
system may be obtained from equation (Bl). The result is 




- it ( c % 6 + c V ) 


(Al) 


For measurement of C^ the gimbal was arranged to allow pitch freedom only. 
Hence, for t = 0, the pitching-moment coefficient is 


^ qSD 


(A2) 


where K^G is the moment measured by the strain gage. The experimental data 
are presented in figures 22 and 23 , where the moment coefficient is plotted 

against the measured pitch angle 0 . The slope of these curves is C ^ . 


Two distances between the gimbal axis and propeller plane were used in 
order to separate the force and moment contributions to the total moment. The 
equations are 




™e,i - Cm Q ~ 


c ®e,2 = c ^e - c 2© 




where the subscripts 1 and 2 refer to data for Z 0 /r = 0.346 and l 0 j\ 
respectively. Equations (A3) are solved to give 


R 


Cm e ~ 1 2 - “ Z 1 C “0,2) 




2R 


l l 


( C %,1 " ^, 2 ) 


(A3) 


0.691, 

(A4) 

(A5) 
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from which C, 


are derived are 


The quantities 
given in figure 
angle. 


and- C 


% 


%,1 and C ne )2 fromwhich % 

11. Figure 10 shows the variation of C^g and 


% 


with blade 


Determination of and 

The static equation governing the yaw deflection of the engine-propeller 
system may also be obtained from equation (Bl) by letting c “e " Cr > and 
* * 

c nty = -Cdq. The result is 

k 2 f = + C^ej (A 6 

For the determination of the gimbal was allowed yaw freedom only. The 

yawing-moment coefficient is found to be 


P ** 


K t t 

qSD 



c ne e 


(A7 


where K^\|r is the moment measured by the strain gage. The experimental data 

/ 

are presented in figures 22 and 23. The slopes of these curves (C^ ) are used 
to obtain by applying equation (A7); that is. 





(a8 


The ratio \|r/9 is a constant for a given blade angle and dynamic pressure, and 
was obtained by plotting t (as determined from strain-gage output) against 0 
(as determined from the sting angle of attack). Since \|r/9 depends on dynamic 
pressure, C** does also, as shown in figures 22 and 23 . The relation 

was used to obtain C^. 

By using = -C^., the resolution of C^. into its components is carrie 

out in the same manner as that for cJL (see eqs. (A3)). The results are 


p„ _ 

1/7 c* 

J 2 c ”e,i) 

(A9 

wity - 

1 2 - Zpl 1 “e>2 

% 

- 2R - (c* 

1 2 - 

C %i) 

(A10 
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The quantities 1 and g from which and are derived are 

given in figure 11. Figure 10 shows the variation of and with blade 

angle. 
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APPENDIX B 


EQUATIONS OF MOTION 


The equations of motion for a two-degree-of -freedom engine -propeller system 
disturbed from equilibrium are derived in reference 2. With some changes in 
notation, the equations for a system symmetrical in pitch and yaw are (omitting 
the aerodynamic acceleration terms) 


e " + !* ,+g ! 0,+k20ss K ( C “e 9 + C “q 9 ' + C V + C “r* ')] 
r - | 0 ' + g | r + + C^t' - C^e - Cm r e ')J 


(B1 


The propeller aerodynamic derivatives are 

Total pitching -moment derivative due to pitch angle: 

^ - S) - k (B£ 

Total pitching-moment derivative due to yaw angle: 

^ = W ~ 2R % 

Total pitching-moment derivative due to pitch rate: 

Cm q = C m q - ^ | (Bi 

Total pitching-moment derivative due to yaw rate: 

Cm r = C mr - ^ C Zr - ^ Cjh^ ( B " 

* 

For the present purpose, C z and C mr , which are considered small in this cas< 
will be neglected. The derivative C^, which was neglected in reference 2, is 
included here. 


The assumption of a harmonic time dependence leads to the characteristic 
equation for neutral stability. The real part yields the whirl frequency 
equation 
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A 2 = P - ^P 2 - Q 


(B6) 


where 


P = 1 + 


2 vji " KCm e 


(t) 


Q = 1 - 2k Cj 


™e(k) 


(B7) 




Terms involving products of aerodynamic and damping terms are small and are 
therefore neglected. The imaginary part of the characteristic equation gives 
the viscous damping required for neutral stability: 


« -‘C 




q 1 \ 2 r* /lX2 

1 - A - kCtttq^ 


m 


(B8) 


The structural damping required for stability is (from ref. 2) 

g = 2£A (B9) 

An approximate expression for the frequency is obtained by assuming 

HI. JxP <<: x 

J k 


The result is (from ref. 2) 


A = 1 


H 

2Jk 


(BIO) 


Since k is proportional to p, and KC me(j“) is small in comparison with 

Q 

1 - A , equation (b 8) indicates that 2£/p should be approximately independent 
of the air density. This result was used to adjust the experimental whirl data 
for the small variations in density. The values of damping coefficient (g or 
2£ ) plotted were obtained by multiplying the average measured coefficients by 

Po/ P ' 
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TABLE I.- SUMMARY OF MODEL CHARACTERISTICS 


I x , slug-ft 2 

Iy (for Z 0 /R = 0.346), slug-ft 2 
Iy (for Z 0 /R = 0.691), slug-ft 2 
> 1 1 

• 75R> n 

a 


O.OO858 

0.0634 

0.0937 

0.8438 

O.1823 

0.1834 
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TABLE II.- MODEL TEST DATA 


(a) Sting-mounted configurations with Z 0 /R = 0.346 


P0.75R' 

deg 

f e> 

cps 

fty , 

cps 

CD 

C\J 


P> 

slugs/cu ft 

n, 

rps 

Vf, 

ft/sec 

f f> 

cps 

35 

9 . 

20 

9.12 

0.0060 

0.0090 

0.00211 

38.4 

120 

6.88 

46 







X 

.00210 

26.2 

118 

7.37 

52 







.0163 

I 

22.0 

122 

7.63 

58 

*> 

f 

> 

r 

\ 


.0126 

1 

19.4 

139 

7.73 

25 

9.71 

9.65 

.0058 

.0131 

.00210 

60.0 

125 

6.20 

46 







X 

I 

28.8 

131 

7.80 

52 







.0076 

X 

24.1 

138 

8.10 

58 

> 

r 

> 

r 

\ 

f 

.0098 

.00211 

20.0 

144 

8.20 

25 

10.43 

10.26 

.0066 

.0216 

i 

64.4 

136 

6.42 

35 








.00210 

40.2 

128 

7.60 

46 





.0060 



I 

30.4 

136 

8.16 

46 





i 



i 

31.8 

143 

8.05 

52 


r 

1 

r 

.0066 

> 

f 

.00209 

26.1 

148 

8.39 

35 

10.96 

IO.96 

.0150 

.0238 

.00217 

51.8 

159 

7.58 

46 









.00218 

37.3 

167 

8.30 

52 









.00209 

33-0 

186 

8.63 

58 

> 

r 

> 

r 

\ 

* 

> 


i 

29.2 

210 

8.75 

35 

11. 

52 

11.50 

.0152 

.0218 

.00218 

56.2 

173 

8.08 

46 









.00217 

41.0 

188 

8.83 

52 










35 -^ 

202 

8.98 

58 

> 

r 

> 



r 

> 

r 

.00216 

30.1 

222 

9.27 

25 

10.44 

10.40 

.0153 

.0263 

.00218 

72.1 

152 

6.26 

33 







i 

1 

53*4 

150 

6.96 

46 





.0126 

.0200 

'i 

36.0 

161 

8.13 

52 





X 



.00219 

29.4 

165 

8.55 

58 

\ 

r 

> 

r 

.0150 

> 


i 

24.2 

171 

8.66 

25 

10.46 

10.51 

.0175 

.0398 

.00219 

81.3 

172 

6.00 

33 







4 r 

.00220 

58.7 

163 

6.62 

46 







.0316 

4. 

41.4 

182 

7.61 

52 







.0398 

.00219 

35-3 

196 

7.71 

58 

N 

f 

> 

r 

> 

f 

.0368 

i 

29.0 

210 

8.19 

25 

10.71 

10.28 

.0371 

.0393 

.00220 

85.6 

179 

5.56 

46 





.0300 



1 

44.6 

196 

7.60 

52 





.0371 



1 

38.4 

213 

7.73 

58 

> 

r 

> 

f 

.0300 

> 


t 

33-2 

231 

8.10 
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TABLE II.- MODEL TEST DATA - Concluded 


(b) Sting-mounted configurations with Z 0 /R = O.69I 


P0.75R' 

deg 

cps 

1 

V 

cps 

2 £ 

0 

2 £ 


p> 

slug/cu ft 

n, 

rps 

— 

Vf, 

ft/sec 

f f> 

cps 

25 

7.60 

7.60 

O.OO56 

0.0053 

0.00226 

59-4 

120 

4.94 

35 









.00227 

48.8 

145 

5.27 

46 









.00226 

37-6 

171 

5-57 

52 


r 

\l 

1 


/ 

> 

V 



31.0 

172 

5.66 

27 | 

7.64 

7.68 

.0060 

.0074 

1 

s 

f 

63.2 

144 

4.73 

35 









.00220 

44.4 

134 

5-50 

35 







.0062 

4- 

46.8 

141 

5-39 

39! 








.00226 

40.3 

140 

5.78 

46 







.0074 

.00220 

33-2 

147 

5.81 

52 







.0062 

. 00221 

28.4 

156 

6.00 

58 

\ 

< 

V 

f 

S 

V 


.00220 

26.3 

183 

6.09 


(c) Wing-mounted configurations with Z 0 /R = 0.346 


P0.75R' 

deg 

f e> 

cps 

cps 

2£ 

0 


* 

p> 

slug/cu ft 

n, 

rps 

Vf, 

ft/sec 

cps 

25 

9.21 

9.34 

0.0072 

0.0132 

0.00232 

61.5 

125 

5.70 

35 









.00231 

41.0 

121 

6.58 

46 









. 00232 

29.0 

126 

7.24 

52 










25.0 

135 

7.40 

58 

> 

f 


S 

\ 

f 

V 


. 00231 

1 

21.5 

148 

7.67 

35 

11.6 

11.7 

.0290 

.0324 

1 

63.5 

188 

7.60 

46 









.00230 

46.5 

202 

8 . 4 i 

52 









.00228 

4 o.o 

216 

8.60 

58 


/ 


f 

si 

f 

s 

/ 

.00233 

36.0 

251 

8.94 
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L- 61-6611 

Figure 1.- Photograph of sting-mounted propeller whirl model in the test section of the Langley transonic dynamics tunnel. 





L-61-291 

Figure 3*~ Photograph of wing— mounted propeller whirl model in the test section of the Langley transonic dynamics tunnel. 




Blade width 
Tip width 



Figure 4.- Distribution of propeller twist angle and blade width for the airplane propeller of this 

analysis. 
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ro 

vo 


Figure 6.- Photograph of sting-mounted propeller whirl model showing details of design. L-6l-66l4.1 


Simulated engine mass (shaker) 


•Gimbal restraint spring 
(fixed to shaft housing on 
center line of shaft axis) 



(a) Power plant, gimbal, and mounting for whirl-flutter model. 
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Figure 9.- Sketch of cantilever-wing model in tunnel, showing measured node lines and frequencies 
of the wing. All dimensions are in inches unless otherwise specified. 
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Experiment 
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Figure 10.- Concluded. 
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Figure 12.- Comparison of experimental and theoretical lift-lag phase angle 
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(e) Po.75R ~ * 

Figure 14.- Concluded 
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Figure 17.- Comparison of experimental whirl-flutter boundaries with theoretical boundaries based 
on various aerodynamic derivatives with viscous damping. Z 0 /R = 0*691; = °*0059* 








Figure 19.- Comparison of C mq inferred from measured whirl speeds with Cm q derived from theory 








PQ.75R > de 9 


Figure 21.- Comparison of experimental whirl -flutter boundaries for the combination of propeller 
and power plant mounted on the sting and mounted on the cantilever wing. lo/ R ■ 0.^46. 
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(c) Yawing moment. Z 0 /R = 0.691 


Figure 22.- Continued. 











